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Grover's algorithm provides a quadratic speed-up over classical algorithms for unstructured 
database or library searches. This paper examines the robustness of Grover's search algorithm to a 
random phase error in the oracle and analyzes the complexity of the search process as a function of 
the scaling of the oracle error with database or library size. Both the discrete- and continuous-time 
implementations of the search algorithm are investigated. It is shown that unless the oracle phase 
error scales as 0(N~ 1 ^ 4 ), neither the discrete- nor the continuous-time implementation of Grover's 
algorithm is scalably robust to this error in the absence of error correction. 



I. INTRODUCTION 



Issues of fault tolerance and error correction are of both theoretical and practical interest in designing and imple- 
menting quantum algorithms. One approach to diminishing the effects of error focuses on error correcting codes, which 
can be used to correct errors once they have occurred Another approach which has also proven successful is 

that of decoherence-free subspaces, which involves mapping a computation onto a subspace that is relatively free from 
error 3j. This second approach is of interest because it is an example of passive, rather than active, error correction. 
The success of the passive approach leads to the question of whether existing quantum algorithms are inherently 
robust to errors, and, if they are not, whether there might exist modified implementations of these algorithms that 
are robust to certain errors, i.e. without requiring active error correction. 

In this paper, we analyze the robustness of Grover's algorithm to error. Grover's search algorithm is the basis 
for a number of quantum algorithms 0, 0, 0. Most of these algorithms involve the use of an oracle, a black box 
device which takes as input a quantum state in the computational basis and returns as output some function of that 
state. Given a perfect, noise-less oracle, Grover's search algorithm attains a quadratic speedup over classical search 
algorithms. However, this speed-up is predicated upon the perfect implementation of the oracle. Although, for the 
purposes of analysis, the oracle is often treated as a "black box" whose inner workings are unknown, any physical 
implementation of Grover's algorithm must also include a physical implementation of the oracle and this may be 
imperfect. Therefore it is of interest to ascertain what effect the accuracy of the oracle implementation has on the 
overall complexity of the algorithm. We make use of this term here as it is employed in computer science terminology, 
namely, the complexity of an algorithm is defined as the number of computational steps required in order to achieve 
a pre-determined fixed probability of success. 

Several previous papers have studied the effects of oracle noise on Grover's algorithm, using various models 0)H, 
Ifj Hfll ITU . We consider here the random phase error model, addressing it within both discrete- and continuous-time 
implementations of the search algorithm. The effect of random phase errors on the discrete-time Grover algorithm 
was already studied numerically in 0- In this paper we derive analytic results for this model and present numerical 
evidence verifying the validity of these results. We analyze the complexity of the search algorithm as a function of the 
scaling of the errors, and arrive at bounds on the error that must be satisfied for a constant probability of success. 
In particular, we find that in order to achieve a constant success probability independent of the library size that is 
being searched, the oracle error must scale at most as iV -1 / 4 where N is the library size. We also analyze the effect of 
phase errors on a continuous-time analogue of the search algorithm. This is found to be relatively straightforward to 
study within a density matrix approach. We find that this leads to similar results for the complexity as a function of 
the scaling of the errors, namely that there is an upper bound of 0{N 1 ^ 2 ) on algorithmic complexity for 5 > 1/4 and 
0(N 1 ~ 25 ) for S < 1/4, when the error scaling goes as N s . We further show that these bounds are tight; in other 
words, that the algorithmic complexity has a lower bound of f^TV 1 / 2 ) for 5 > 1/4 and ^(TV 1-25 ) for S < 1/4. (Following 
standard computer science notation, we will use the following to characterize the growth of certain functions: We will 
say that f(n) is bounded above by g(n), i.e. that f(n) — 0(g(n)), if there are positive constants c and k such that 
< f(n) < cg(n) for n > k. Similarly f{n) is bounded below by g(n), i.e. f(n) = Q(g(n)), if < cg(n) < f(n) for 
constants c, k > and n > k.) 
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This result has practical importance in determining the amount of oracle error that is allowable as the size of the 
library is increased. We will show that our complexity results imply that if the library size is increased by a factor of 
k, the oracle error must decrease by a factor of fc 1 / 4 in order to attain a similar probability of success. 

The remainder of this paper is organized as follows. Sec. ITT1 describes the random phase error model for a noisy 
oracle, summarizes the Grover search algorithm in a discrete-time implementation, and then derives the effect of the 
oracle noise on this implementation. Sec. IIIII derives the effect of oracle noise on the continuous-time formulation of 
Grover's algorithm proposed by Farhi et al. in p"2lll3j|. We summarize and present conclusions in Sec. II VI 



II. DISCRETE-TIME SEARCH 

A. Quantum Search with a Phase Oracle 

A phase oracle is a quantum oracle that "marks" one or more computational basis states with a specified phase 
(usually —1). For a function /(x) — > {0, 1}, where x denotes an n-bit binary string, a phase oracle implements the 
following operation: 

|x)^(-l)/( x )|x) (1) 

The search problem is phrased as follows: Given an oracle, O, which implements the function /(x) — * {0, 1}, find a 
state in the computational basis x for which /(x) = 1. Assuming a perfect, noiseless oracle (the concept of a noisy 
oracle is introduced below), the well-known result for the complexity of Grover's algorithm, i.e., for the quantum 

search problem, is O(y^j^), where M is the number of states for which /(x) = 1 [T3. ITU llil Il7 | . In contrast, the 

classical complexity of the search problem is O(-p0 |18| . Thus, the quantum algorithm provides a quadratic speed-up 
over the classical algorithm. 



B. Grover's algorithm in absence of noise 

The description of Grover's algorithm in this paper follows the discussion in The search is made on the set 
of N = 2™ states represented by all n-bit binary strings x £ {0,1}™. The states are denoted by \x) . Within a 
discrete-time computation model, Grover's algorithm searches for marked or "target" states via repeated application 
of the Grover iteration operator, G, which can be written as: 

G=(2\r,)(r,\-I)0 (2) 

where = X^^Lo 1 l x ) i s the equal superposition over all states. From this point on, we will assume that M = 1 
to simplify the discussion. The arguments can easily be extended to the case where M > 1. Let the state |r) be the 
"marked" state. Then we can write out the explicit form of the oracle as: 

= 7+(e--l)|r)(r| (3) 



Assuming that our registers are initialized to it can be shown that after 0(vN) applications of G, the quantum 
computer is approximately in the marked state |t) 0,0]. This result is demonstrated by noting that the search 
algorithm can be mapped onto a two-state subspace spanned by the two basis vectors |1) and |2), where 

|1) = |r>, (4) 
|2) = (5) 



In other words, |1) is the marked state and |2) is the equal superposition over all unmarked states. 
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When our initial state, \rj), and the operator G are rewritten in the |1), |2) basis, 

\V) = [ y^rl (6) 




r - 1 " W \ ( cos(e/2) sm(6/2) \ 

where 

sin(6/2) = a^EI, (8) 

we see that the effect of repeated applications of G is to make successive rotations by 9/2 on state |?7). This convention 
for the definition of 9 is useful because 9 represents the angle of rotation applied by the operator G on the Bloch 
sphere defined for the Grover subspace spanned by the basis states |1) and |2). We can therefore view Grover's 
algorithm as the process of rotating our initial state, |fy), onto the target state, |r) = |1) by a discrete sequence of 
small rotations. For small 9 (i.e. large N), sin 9 « 9, yielding an incremental rotation angle of approximately 
^n^ 1 ~ Then a rotation of n radians on the Bloch sphere requires O(VN) applications of G. Hence, Grover's 



search algorithm requires 0(v TV) calls to the oracle. The well-known lower bound for quantum search has been 
established as tt(V~N) [2(|. Thus, Grover's algorithm is an optimal search. 



C. Grover's algorithm in presence of oracle noise 



This 0(^/N) bound on the running time of the search assumes that the oracle exactly implements the phase 
operation specified by Eq. Instead, we could envision a more realistic oracle which implements Eq. with some 
finite degree of precision. Specifically, we adopt here the concept of a noisy phase oracle introduced in |]j which 
implements the following operation: 




where e is a small, random phase error (e <C it) with probability distribution p(e) |23j |. We will make the assumption 
that the error is zero-mean, i.e. that (e) = 0. For p{e) — (5(e), the noiseless oracle is recovered. In practice, the root 
mean-square magnitude of e can be made small through a careful physical implementation of the oracle. However, 
due to the finite precision of any experimental implementation, the average magnitude of e will never be zero, i.e., 

Given this more realistic oracle model, we now investigate whether Grover's search algorithm is still 0{-jN). The 
explicit form of the noisy oracle can be written as, 



= /+(e^ +e >-l)|r) (t\ 
Then Grover's iteration operator, G, becomes: 



G 



(-1 + AW(w+e) zV^T 

2yOV-T i(n+e) i _ 2_ 
N c N 



(10) 



(11) 



which is clearly still unitary. To sec the effects of the random error, we first note that any single qubit unitary operator 
can be written as 



U = cxp(ia)Rn(&) — exp(ia) cos — ] I — isin I — In • a I , 



4> 



(12) 
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where a are the Pauli operators. Thus, the Grover's iteration operator using the noisy oracle is equivalent, up to an 
overall global phase factor, to a rotation, i?^(<I>), on the Bloch sphere by some angle $ about some direction h. We 
can determine the value of $ and n by using Eq. (flip and Eq. (|12|l . 

$ 

n,S| 2 



$^ 2\/N — 1 

sin 



«MI) < i3 > 
(I) <»> 



2/ TV 

sm^-j % = — (15) 

rin(|)n, = -(l-f)sin(f). (16) 

A noiseless oracle can be recovered using these equations by setting e to zero. In that case we obtain the large N 
Grover rotation angle, $ w with n x = rt z = 0, and n y = 1. Hence, the noiseless Grover's algorithm corresponds 
to a rotation on the Bloch sphere about the y-axis. The effect of a noisy oracle is to add small, random perturbations 
to this rotation axis, thereby changing the pure y-axis rotation to a rotation about an axis n which also contains 
non-zero x- and z- components. 

To analyze the running time of the noisy search algorithm, we will take the limit of large N. This is a useful 
assumption for our current purposes since we are interested in the complexity of the search algorithm for large N. 
From Eqs. (|13H16|) . we can express G up to a global phase factor as, 



G = exp ( - « ( ( -1 + |) sin (|) a z + (sin (|) a x - cos (f ) a y ) ) / sin (f ) ) , (17) 



where sin(^) = y 1 — (1 — -^) 2 cos 2 (|). Since both e and N 1 I 2 are small parameters, we can use a double Taylor 

expansion of the terms in the exponent with respect to e and iV -1 / 2 . Excluding second-order terms and higher, we 
obtain 

(■( 

G = cxd I i I eov H ; 

N 

Again using the fact that e and TV -1 / 2 are both small parameters, we can make use of the Baker-Campbell-Hausdorff 
formula [23 to arrive at 



G = exp ( i ( ecr z H — 7=0^ ) ) • (18) 



G = exp (iea z ) exp \i^=a y ^j + 0{^) 
« J?i(-e)i^(-^). 



(19) 



Examining Eq. I|19|l . we observe that there are two important timescales that will be relevant to any physical 
implementation. The first is the timescale of the y-axis rotation, y/~N. Taking the limit as e — > 0, we see that the 
y-axis rotation, R y is the "driving force" of the search algorithm, in that it rotates the initial state, |?y), towards the 
target state, |r). The second important timescale, 1/e, is the timescale of the random rotation about the z-axis of 
the Bloch sphere, This rotation is the effect of noise and does not advance the search process. Because e is a 
random variable, this rotation about z will in general be different for each application of G. However, e does have a 
characteristic root-mean-square magnitude, e rms , which is determined by the probability distribution p(e). Then the 
characteristic timescale of the z-axis rotation is determined by l/e rms . 

In general, e rms and ^/N are independent parameters: e rms is the phase error that is specific to the experimental 
implementation, while N is the size of the database. However, in reality, these two parameters will be linked for a 
given experimental setup because the size of the database (i.e., the number of states N — 2 n accessed by the n qubits 
of the computer) will affect the accuracy of the oracle implementation. Note that the oracle acts on all n qubits, (see 
Eq. ||SJ|). It seems very unlikely that e rms would decrease as N increases, since a greater number of qubits generally 
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introduces more potential for noise. In the best case scenario, it might be possible to find a physical implementation 
for which e rms is essentially constant over a large range of values for N . To ensure that our analysis is as general 
as possible, we will assume here that e rms scales as N~ s , where S is some constant that may take both positive and 
negative values, or zero. We can then determine what effect the scaling constant S has on algorithm performance. In 
particular, we shall determine the effect of S upon the complexity of the algorithm, i.e., upon the number of Grover 
iterations required to achieve a given probability of successful measurement of the target state This will allow 
us to further determine the maximum allowable oracle error scaling that ensures a constant probability of success 
independent of the library size, for a given rate of library growth. 



D. Dynamics and scaling of the noisy Grover search algorithm 



To analyze the dynamics of the Grover iteration, we observe that any pure state on the Bloch sphere can be 
represented in spherical coordinates by two angles, 9 and 0, where we define 9 to be the angle from the z-axis and cf> to 
be the azimuthal angle. Here we have chosen our coordinates such that the target state |1) is the south pole (8 = it) 
and the state |2) is the north pole (9 = 0) of our Bloch sphere. Because the operator G depends upon the random 
variable e, each application of G induces a new probability distribution over the Bloch sphere which we will denote 
by f{9, <p). In this notation, our initial state at time t = is given by a delta function (i.e. a pure state) centered on 
= 0, with 9 = sin _1 (l/\/]V) w 1/VN. 

Using this notation, we can write the probability of obtaining the state |1) upon measurement after t iterations as 

P(t) = 1- ± J Jcos 9 f t (9,(j)) sin 9d9d(/) 

= i(l-<cos0» (20) 

= i(i-<*». 

Let us define the complexity of our search algorithm to be the number of iterations required to obtain a success 
probability of 2/3. We immediately see that if the probability of obtaining state |1) upon measurement after t 
iterations is P(t), then by repeating this procedure approximately gp^y times, we can boost the overall success 
probability to 2/3. It should be noted that the choice of the constant 2/3 is arbitrary; in general, the choice of 
constant will not affect the complexity of the algorithm. Thus, the overall complexity of our algorithm is 0(t/P{t)). 
Furthermore, Eq. (|20|l states that in order to succeed with some desired probability P*, our probability distribution 
function f{9, <f>) must be non-negligible when 9 > 9* , where 

9* = cos" 1 (I - 2P*) . (21) 

In other words, to obtain some desired probability of success, P* , there must be a high probability of reaching points 
on the Bloch sphere with polar angle greater than 9*. 

Having established this terminology, we will now give a phenomenological description of the evolution of the prob- 
ability distribution over the Bloch sphere. Specifically, we would like to know what dependence the magnitude and 
scaling of the error e rms has on the maximum attainable polar angle. 

We first rewrite the effects of the Grover iterate as a function of polar coordinates. In the polar coordinates defined 
above, the z-rotation can be written as 

<f> 5$ + e , (22) 
9 5$ e. (23) 

If we momentarily neglect the y-axis rotations, our z-axis rotation dynamics correspond to a random walk on the 
variable <f> with periodic boundary conditions. To analyze the effects of the y-rotation, we then take advantage of the 
fact that 1/y/N is a small quantity and expand in powers of 1/yN to obtain 



cos^+Ofl/W). (25) 
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Finally, we can write the combined effects of our noisy Grover operator G using Eqs. I|22H25|I . as 

, g , . ,cos# 4 

cj> -> 0-sin0-— -= + e, 26) 
smfc> y^/V 

6> 6/ + cos^^, (27) 



where we have dropped terms of 0(1/ N). 

Having written the dynamics of the Grover operator in terms of polar coordinates on the Bloch sphere, we now 
consider the probability distribution of the quantum state over the Bloch sphere after T Grover iterations. We 
will analyze the dynamics for <fi -C 1, and then consider when this approximation is valid. In this regime, we can 
approximate the Grover operator dynamics as 

<f> $ + e , (28) 
N 



(29) 

It is evident that in this small <j> limit, the </> dynamics are completely determined by z-axis rotation, which results in 
a random walk on the variable <f> that is controlled by the random variable e. The central limit theorem tells us that 
regardless of the probability distribution p(e) from which the random variable e is sampled, after an adequate number 
of iterations the distribution of will converge to a Gaussian with width e rms \T: 



,y/T, (30) 



where T is the number of iterations. 

Let us now consider when this approximation is valid. Clearly, the condition that (j) <C 1 is satisfied for the initial 
state (po = 0. We recall that e rms is assumed to scale as N~ s , where S is some constant. The validity of our small 
(j> approximation is found to be highly dependent on the scaling exponent S. We consider two cases: 5 > 1/4 and 
5 < 1/4. 

If 8 > 1/4, we examine the probability distribution after T = X\/N iterations, where A is some small constant. 
Inspection of Eq. (|30|l shows that the probability distribution of <j> after T steps will then have width y/XN 1 ^^ 6 . Since 
iV can certainly be made arbitrarily large, this justifies our assumption that <f> is small. Within this regime, Eq. 128(1 
shows that the dynamics of the variable 9 are simply those of the deterministic y-rotation with constant increment 



-$=■ After T iterations of this, we obtain 

VN 



9 ° + S' 

4T 



(31) 



which is the same as the large N limit of the noiseless Grover search, Eq. |JS}. Thus when 8 > 1/4, the error has no 
effect on algorithmic complexity for significantly large N, yielding a complexity of 0(N 1 / 2 ). This bound is trivially 
tight since the search problem is well-known to have a lower bound of f^TV 1 / 2 )!^. For the remainder of this section, 
we will therefore analyze the complementary case, i.e., when e rms scales as N~ s with 6 < 1/4, and determine the 
effect of the error on algorithmic complexity. We note that this latter case includes the best-case physical situation 
of e r rns independent of N, i.e., 6 = 0. 

In order to analyze the system dynamics for S < 1/4, let us examine the probability distribution after T iterations 
such that 

T = A/ e 2 rms , (32) 

where A is a small constant. By Eq. 1(30(1 . the distribution width 4> rms is proportional to vA. Thus, by selecting a 
small enough constant A, our approximation that (j> is small is again justified. Using Eq. 128(1 and Eq. 1(31(1 . we find 
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that the ^-rotation is again essentially deterministic and that 

4T 



N 
4A 



rms 



Using Eq. ((201), we obtain 



N 

= (,(—L-\. (:-!:-!) 



(34) 



where we have used the fact that - x ^ <C 1 when S < 1/4. By proving that we can attain (essentially de- 
terministically, as described above) a polar rotation of at least - 4- ^jy in T = \/e 2 . ms iterations, we have shown 
that the complexity of the noisy search algorithm is O(p^fy) — 0(e 2 . ms N). Taking the scaling of e rms into account 

(e rms ~ N~ 5 ), leads to the overall algorithmic complexity 0(A^ 1 ^ 2 ' 5 ) 24]. So for 5=1/4 we obtain an optimal solution 
having the same speed-up as the noiseless quantum search. For < 5 < 1/4, we find a scaling intermediate between 
the noiseless quantum and classical search algorithms. The quantum speed-up factor decreases as S approaches zero 
and is completely lost when 5 — 0. (Note that the quantum search formally becomes slower than the classical search 
in the worst case scenario when the phase errors increase with N, i.e., S < 0.) 

In order to show that this complexity bound is tight, we now will show that we can attain at most a polar 

rotation of O ^ - ^ in sub-classical time. Repeating the argument above following Eq. (|33[l will then lead to the 

identification of ^iN 1 - 25 ) as a lower bound. Thus, we find that 0(N 1 2S ) is a tight bound on the complexity in 
presence of noise. Obtaining a lower bound on the algorithmic complexity is important because without it, it would be 
unclear whether we could further reduce our search complexity by running the noisy search algorithm for more than 
0(1/ e^ms) iterations. The lower bound will demonstrate that after 0(l/e^. ms ) iterations, it is impossible to achieve a 
super-classical enhancement in success probability. 

To achieve a polar rotation greater than O ^— — ~r^j j we will need to obtain a significant probability density on the 
region of the Bloch sphere described by 9 > e2 c ^ , where c is a constant. Let us therefore consider the dynamics 
of a state on this region of the Bloch sphere. Given a pure state with (f>o = and ^ > 9 > - we examine the 

action of T iterations of the Grover algorithm on this state, where T = \/t 2 ms . Now, instead of choosing A small 
as previously, we select A > 1, which ensures a large width for the long time distribution of the azimuthal angle <fi, 
Eq. (|30(l . Thus after T iterations with this large width, we obtain a completely uniform distribution on the variable 
(\>. Wc now consider what happens to the polar variable, 6. Referring to Eq. (|23|l and Eq. I|25|l . we find that after T 
iterations we obtain 



\e T -e \ < 1 



< 



O (7V- 1 /2+25) 



(35) 



Again, provided that 8 < 1/4, this difference can be made arbitrarily small by selecting a large enough N. Thus, 
after T iterations, our probability distribution of the polar coordinate 9 is confined within an arbitrarily small window 
around 9o, but is uniformly distributed on the azimuthal coordinate 4>. Figure 1 shows a graphical representation of 
these dynamics. 

Let us now consider the effects when the Grover operator is subsequently applied to such a probability distribution 
that is uniform in <f>. We will show that there exists a symmetry in the associated 9 transformations in this regime, 
such that one half of all azimuthal angles, <f>, are associated with a 9 rotation in one direction, and the other half are 
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associated with 9 rotations of equal magnitude in the opposite direction. To reveal this symmetry, we show that the 
transformation cj> — * (f> + 7r on Eqs. (|26|) and (|27|) results in an exactly opposite rotation in the 9 direction to that 
associated with the azimuthal angle, (p. From Eq. (|35|1 . we know that the distribution of 9t will be arbitrarily close 
to #o- Since we have assumed that our initial state was located in the region tt/2 > 9q > £-^=, we can rewrite Eq. 

(|26|l as follows: 

, o . cos 9 1 



sin6» y/N 



■2 



N 1 

+ € 



C y/N 

,_ sm( p™i< +e . ( 36 ) 



We now make use of the arbitrariness in choice of the constant c, choosing c 3> 1 so that c 3> sin 0, making the second 
term arbitrarily small. Hence we obtain 

4> $ + (37) 

which demonstrates that the transformation <j) — > <p + 7r has no net effect on the azimuthal dynamics. On the other 
hand, our 9 dynamics become 

9 ^ 9 + cosU + n)^ + 0(l/N) 
WN 



6>-cos0— L + 0(l/iV). (38) 

Wn 



Eq. (|38(l shows that the transformation <fi — > <f> + ir causes the polar angle 9 to be rotated in the opposite direction, as a 
consequence of the change in sign of the cos0 factor. It follows that once cf> has reached a uniform distribution, there 
will be an equal probability of rotating 9 by A/^/N in the positive and negative directions. This behavior is valid on 
the entire region of the Bloch sphere specified by our initial condition tt/2 > 9 > - r -^=. As a result, the 9 dynamics 

over this entire region no longer resemble a deterministic rotation towards 9 — tt, but instead resemble a random walk 
on 9 with step size of approximately A/\/N. Figure 1 shows a graphical representation of these dynamics. 
Using the central limit theorem, we find that 

9 rms cx y/t/N, (39) 

where t is the number of iterations after the uniform distribution in <fi has been attained, i.e., after the first T 
iterations. Consequently, to "move" the polar angle of the probability distribution by some small angle A9 now takes 
time proportional to A9 2 N. From Eq. (J2DJ, the corresponding change in the probability of a successful measurement 
of the target state is 

AP w - sin A9A9 « i A9 2 cx At/N (40) 

Thus, once we have entered the region of the Bloch sphere characterized by tt/2 > 9 > - the probability 

of a successful measurement further increases only at the rate t/N. This is the same as the classical result and 
consequently, we can attain no further speed-up over the classical algorithm once we have entered this regime. The 
polar angle 0( ^ 2 ^ ) thus constitutes an upper bound on the rotation that can be achieved in sub-classical time. 

Hence, by the arguments given above, our bound 0{N 1 ^ 2S ) on the complexity of the noisy quantum search algorithm 
is a tight bound. 
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E. Summary of discrete time noise Grover search 

Using the above results, we can now give an accurate phenomenological description of the dynamics of the noisy 
search algorithm. We identify the critical timescale of the algorithm as t m i X i ng = l/e 2 ms , which is the mixing 
time of the z-rotation induced by the random error. For times t <C t m i X i ng , the algorithm proceeds more or less 
deterministically, and the initial state is rotated towards the target state |1), attaining a polar angle of 8 = e2 1 v / ]y • 

However, for times t 3> t m i X i ng , the <j) variable becomes completely mixed and the ^-rotation, which is the driving 
force of the search process, has an equal probability of increasing or decreasing the polar angle of the probability 
distribution. As a result, the 6 dynamics also become those of a random walk, and the searching rate becomes 
classical. Therefore, in order to gain the maximum speed-up over the classical algorithm, the best methodology we 
can employ is to run the algorithm for t <C t m i X i ng , measure, and then repeat this process to boost our overall success 
probability. This protocol will allow the sub-classical scaling to be retained. 

The effects of the scaling of e rms on the complexity of the search algorithm are seen to be different according to 
the value of the error scaling index 5. We saw that if S > 1/4, the search algorithm was essentially unaffected by the 
presence of error and the complexity is identical to the noiseless quantum search result of 0(y/N). On the other hand, 
if 6 < 1/4, then the optimal speed-up is obtained by running the algorithm for t = 0(l/e 2 ms ) iterations, achieving a 
maximum polar angle of 9 = 0(l/e 2 ms \/N) and then measuring. Expanding Eq. (|20[1 yields the measurement success 
rate, P{t) = 0(l/ef. ms N) and a resulting algorithmic complexity of 0(t/P(t)) = 0(e 2 rm8 N) = OiN 1 ' 26 ). So the 
optimal error scaling is given by 5 > 1/4, i.e., the errors scale as 0{N~ 1 ^). In contrast, the most physically realistic 
constant error scaling 5 = 0, corresponding to a constant error over a range of N values, yields no speed up over the 
classical search. 

It is interesting to compare these results to the analytic results for a constant phase error oracle given in • In the 
constant phase error model, the oracle applies the same phase exp + e)) to the marked state at each iteration. 
It can be shown that for a constant phase error of magnitude e, the error must scale as £ = 0(N~ 1 ^ 2 ) in order to 
obtain a quadratic speed-up over the classical search. The differences in error accumulation between the constant 
and random phase error processes can be compared (for example) to the ballistic and diffusive regimes of Brownian 
motion, respectively. Constant phase errors tend to accumulate quickly since they are all in the same "direction" . On 
the other hand, subsequent random phase errors can cancel each other out and thus accumulate more slowly. 

Let us now consider a question of physical importance. Given an oracle with error magnitude e rms and a library 
of size N, let us assume that we can attain some constant success rate P after T iterations. Given a larger library 
of size N' = kN, we wish to obtain the same success rate P after T' = y/kT iterations. In other words, we wish to 

obtain the ideal quantum search complexity of O ^y/W^j also for a larger library size. What is the maximum allowable 

oracle error, e' rms , for the larger library? The answer to this question follows immediately from our previous analysis. 
We know that e rms can scale as N~ s with 6 > 1/4 without affecting algorithmic complexity. Thus the maximum 
allowable error in the oracle for the larger library is e' rms — fc _1 / 4 e rms . As an example, if the library size is doubled, 
the oracle error for the larger library is required to be at most 0.84 the corresponding error for the smaller library. 

To verify the validity of our analytic results, we can simulate the effects of a noisy oracle numerically. Given a 
library of size N, we simulated the Grover's search algorithm using the procedure described in Section Hi CI For each 
iteration of the search algorithm, a phase error e was selected from a zero-mean Gaussian distribution with standard 
deviation e rms . The algorithm was run for t = n^/N / 4 iterations and the maximal probability of success attained 
was recorded. This process was repeated 100 times and the average success probability was calculated. Figure 2 plots 
the average probability of success, (P), versus the library size, n = log N, for various values of e rms . Figure 3 shows 
the value of loge rms that yields a constant success probability of P — 1/2 for a given library size. As predicted by 
our analysis above, we find that in order to obtain a constant success rate as a function of n (or N), the error must 
scale as e rms cx TV -1 / 4 . 



III. EFFECT OF NOISE ON CONTINUOUS TIME ANALOGUE OF QUANTUM SEARCH 

In this section, we analyze the effect of noise on a slightly different model for quantum search that has been proposed 
by Farhi and Gutmann 12j. In this model, one again starts in the symmetric superposition of all states, \rj), and then 
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applies the following Hamiltonian 

H Q = \v)(v\ + \r)(r\=H v + H T , (41) 

where |r) is the marked state. Note that this Hamiltonian is directly related to the Grover's iterate, Eqs. J2J an d 
(0. Action of Eq. (J2J) amounts to simply applying the Hamiltonian H n for a time ir, followed by applying the 
Hamiltonian H T for a time 7r. Clearly the two operators Eq. (J2J) and Eq. ( 141ft would be equivalent if H v and i/ T 
commuted. However, \H n ,H T ] = -^(|r?)(r| + |t)(t/|). Thus the two methods become similar for large N and are 
formally equivalent as N approaches infinity. 

Farhi and Gutmann calculated the time evolution of the system when one starts in the equal superposition state 
| rj). They found that 




They noticed that at time t — ^ ^ , the initial state, |ry), has evolved to the marked state, |t) = |1). The time 

required to evolve to the marked state scales as 0(\/~N), matching the complexity of Grover's algorithm with respect 
to an oracle. As a result, we take the time that it takes to reach the marked state as a function of N to be the measure 
of the complexity of a continuous time algorithm. 

As noted above, in the limit of large N the continuous time and discrete Grover's algorithm are formally equivalent. 
Therefore, it seems useful to also evaluate the effect of a fluctuating "oracle" in the continuous time picture. The 
continuous time Hamiltonian that models a discrete quantum search noisy oracle with phase fluctuation e is given by 



H = H V + (1 + 0H T 

= |r?)(^|+(l + 0|r)(r|. 



(43) 



Here £ is a time-dependent, Markovian stochastic variable that satisfies J* £dt — e. We shall assume that e fluctuates 
and that it can be described by a Gaussian distribution. 

In order to evaluate the Hamiltonian and its effect on the initial state, it is simpler to transform to our orthonormal 
two state basis |1) and |2) defined by Eqs. (@J and 10. One can then decompose the transformed Hamiltonian into 
the corresponding spin operators to find that 



H 



:i + - ) t |— + ~)« 



l 



i 



N 



1 c T ■ 

N 



(44) 



When the fluctuations in the time evolution operator e %m are now considered, i.e., the resulting phase e ^ 0, it is 
useful to approach the problem by examining the evolution of the density matrix. Consider an initial density matrix, 



p{0) 



a b 
b* d 



(45) 



It is well established that a fluctuating a z component of the Hamiltonian leads to dephasing. For example, applying 
the fluctuating perturbation |cr z for a time ir yields, 

p(t) = exp (-icr z 7r/2)p(0) exp (ia z ir/2). (46) 

Subsequently averaging over all possible values of the resulting total phase, e, leads to the following evolution: 



1 



2tts 



a b 
b* d 



a be- 3 '' 2 
b*e- s2 / 2 d 



(47) 



Here s 2 is the variance of the fluctuating phase e. Since it has been assumed that the fluctuating field £ is Markovian, 
one can treat this as a dephasing term described by a decay constant T = s 2 /2ir and then use the corresponding Bloch 
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equations [22] to calculate the evolution of the system. This Markovian approximation is valid when the fluctuation 
of £ is much faster the 
Bloch representation, 



of £ is much faster than -^=. To analyze the system evolution, we therefore decompose our density matrix using the 



p = \l + \n ■ a. (48) 

The time evolution is calculated by solving the Bloch equations for the components of n : 

2 

n x = j^n y - Tn x , (49) 



W v ^K~l^- rr ^ (50) 



N V V N 



(51) 

In order to understand the effect of the dephasing on the algorithmic complexity, we examine these equations in 
the limit of large N and keep only terms that are of order -^=. This yields only two coupled equations, in n y and in 



2 N n z - Tn v (52) 



n z = j= n v (53) 



N 



Our initial density matrix is given by 



p = \v)iv\ 



1 T 1 / 1 2 \ 1 (2s/N=\\ 



2 2V iV/ 2 V iV 



We note that the quantity n z is the projection onto the state |1) and thus provides a measure of how well the 

computation is proceeding. Initially we have n z = — 1 + while our computation becomes complete when n z = 1. 

Solving the differential equations in the large N limit, (Eq. J33J)), with initial conditions n y (0) — 0, n z (0) — —1 + jj, 
yields the following solution: 



»vW = SSl ( ex P (l (-^- Vr 2 - 16/jv)) +ex P (| Vr 2 - W^v))) 

= ( (~^™ + ^ N - 16 ) ex P (l - Vr 2 - Wiv) ) (55) 

+ (Vr^/v + vT 2 iv- 16) ex P (l + ^r 2 - w/n^ . 

Ideally, we would like to find the time at which n z (t) — 1 , or equivalcntly, the time when the probability of reaching 
the marked state is unity, 

P{t) = {r\p{t)\r)= l -{l + n z {t)) = l. (56) 

When measuring the complexity of Grover's algorithm however, we need only to find the time required such that the 
probability Pit) of being in the marked state r is greater than some constant. For concreteness, we choose here the 
minimum time satisfying P(t) — (r|p(t)|r) > 1/4. 
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To determine this time explicitly, we calculate P(t) as a function of various values of V in the limit of large N. 
We observe that in Eq. (JSSJ the term ^T 2 - 16 /N is imaginary when T < Consequently, we expect that the 

behavior will be drastically different for the two regimes a) T < -7=, and b) T > -^=. In order to make a direct 
comparison between the continuous time behavior and the discrete time results, we choose the dephasing constant T 
to scale with N in the same way as in Section m namely 

r = aN~ 2S cx e 2 rms , (57) 

where a is a proportionality constant. This provides the contact point of the continuous time search algorithm with 
the discrete time algorithm of Section |n] We then see that regimes a) and b) correspond to the two regimes already 
established in Section [H] i.e., in region a), 6 > 1/4, and in region b), 6 < 1/4. We note that our regions b) and a) 
both include S = 1/4. The boundary line between the regions occurs when (5=1/4 and a = 4. 
For T = aN~ 2S < -4=, i.e., in region a), we calculate from Eq. l(S3|l that 



P(t) = \ + \ [-1 + |) e'/ 2r Ls(v/l6/7V-rV2) + V r^jV-ie sin(Vl6/JV-r^/2)j . (58) 



We now pick a time 



and find that 



Hence, 



f = ; (59) 

^16/iV-r 2 



Pit') = - + - (l - -) exp ( , nT 1 • (60) 
V ; 2 2^ NJ v \ ^/16/N -V 2 V ' 



P(t') > \-0(l/N). (61) 



Since this is larger than the value corresponding to our definition of minimum time, it implies that the complexity of 
the search algorithm is bounded from above by 0(t'). Inspection of Eq. I|59|l shows that t' is an increasing function 
of r. Therefore, in order to evaluate an upper bound for t' in the regime a) where 5 > 1/4, we need to evaluate t 1 
for the largest possible value of T. The largest value of T in this regime lies on the boundary with regime b), namely 
where 5 = 1/4 and a = 4. Hence, we evaluate t' as one asymptotically approaches the boundary between regions a) 
and b), for a given error scaling 5. We choose (5 = 1/4 and a = 4 — m, where m is a small constant greater than zero. 

Substituting in Eq. l|tj(J|) yields t' = 2n vN . Thus the minimum time will be on the order of O(VlV), corresponding 
to an upper bound on the algorithmic complexity of 0(^/N). Hence for 6 > 1/4, in regime a), the continuous time 
search algorithm achieves its maximal algorithmic quantum speed up. This is in agreement with the results from 
Section [HI which showed that the discrete time quantum search algorithm yields the maximum algorithmic speedup 
when the root mean square error e rms is smaller than iV -1 / 4 , corresponding to 5 > 1/4. 

To complete the continuous time analysis, we solve for P(t) in the regime b) where S < 1/4. Here evaluation of n z 
in the limit of large N leads to 

P{t) = \ \e-^'^. (62) 



Setting P(t') = 1/4 and solving for t' yields 

, ATln(2) 



4 

rl-2<5 



(63) 



0{N L - za ). (64) 
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We have verified our conclusions by numerically simulating Eq. (|55|l for various values of N and T. Figure 4 shows 
a logarithmic plot of the minimum time t' required to obtain a success probability of P(t') = 1/4, as a function of 
the error scaling parameter S, for various values of N. It is evident that the algorithmic complexity shows a marked 
transition at 5 = 1/4, from scaling OiN 1 / 2 ) for 5 > 1/4, to 0(N 1 ^ 2S ) for 5 < 1/4, as predicted by Eq. (Si and Eq. 
(10411 respectively. Equivalently, we can state that the continuous time quantum search algorithm with randomized 
phase error achieves minimal complexity and hence maximum algorithmic speedup when e rms < iV 1 / 4 (5 > 1/4). 
This agrees with the results for the discrete-time algorithm. 

IV. CONCLUSIONS 

The analysis in this paper has provided a phenomonological description of the dependence of the algorithmic 
complexity of Grover's algorithm on the scaling of oracle phase error for a discrete quantum search, and on stochastic 
Hamiltonian errors leading to phase error in a continuous time quantum search algorithm. In both the discrete- and 
continuous-time versions of the algorithm, it was found that if the phase error scaled with size as N~ s , then for 
6 < 1/4 the effect on the complexity of the algorithm for large N was negligible. However, if the size scaling of the 
error lies in the regime S < 1/4, then for large N it was determined that there is tight bound of O (iV 1 ^ 2 ' 5 ) on the 
complexity of the search algorithm. In particular, this implies that in the presence of any constant (non-zero) amount 
of phase error in the oracle ((5 = 0), there exists some library size N max above which the quantum search algorithm no 
longer provides a quadratic speed-up. In this case, for databases of size N > N max , the search time is il(N), which is 
equivalent to the classical result and there is therefore no quantum speedup. Intermediate error scaling, < 5 < 1/4, 
provides speedup intermediate between the classical and quantum limits, respectively. These results hold for both the 
discrete-time and continuous-time quantum search algorithms, and assume very little about the specific form of the 
underlying error processes. 

The complexity analysis we have made here is also important for determining the precision needed in scaling up a 
quantum search. For instance, let us assume that we are able to implement a quantum search for a library of size N 
with an oracle error of magnitude e rms . Then to perform a quantum search on a library of size kN with equivalent 
accuracy, our results imply that we need to implement an oracle with an error of at most erms/k 1 / 4 . Since this 
must lie in the regime S > 1/4, physically, this requires a system where the phase error decreases as a function of 
database size. Consequently, the precision must increase exponentially as a function of the number of qubits, putting 
severe demands on the physical realization. In contrast, if the error and hence the precision is constant in the system 
size, (e.g., a system where the natural line width is independent of the number of states), then there always exists a 
database size such that the quantum approach offers no speed up over the classical search algorithm. 

The main consequence of the non-robustness of these forms of quantum searches to oracle noise that was demon- 
strated here is a practical limitation on the size of the library on which may be searched with a quadratic speed-up 
using quantum search algorithms without any explicit error correction. This result has significant consequences for 
physical implementation of quantum search algorithms, since although quantum error correction can be used to reduce 
the error present in the oracle, such error correction procedures can require significant resources 2j. In practice it will 
therefore be necessary to balance the cost of error correction (in both spatial and temporal resources) with the extent 
of speed-up attained by a noisy quantum search. The analytic results presented in this paper provide a useful bound 
on the maximum oracle error permissible if a quadratic speedup is to be retained. Above this maximum allowable 
error, we must use error correction in any physical implementation. Conversely, below this maximum error, we can be 
confident that error correction will not be necessary, provided that oracle phase error is the primary source of error. 

As a final comment we point out that although Grover's algorithm and its continuous time analogue are not 
inherently robust to phase error in the oracle, it is not clear whether other implementations of quantum search may 
be inherently robust. Exploration of both active and passive error correction schemes for Grover's algorithm will 
therefore be a valuable direction for future work. 
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FIG. 1: Graphical representation of the noisy oracle Grover algorithm dynamics when J > 9 > ^^/W' The target state |1) 

is located at the south pole of the Bloch sphere. The north pole is state |2), Eq. JjJ. After T iterations, a pure state at 9o and 
<f>o will become completely mixed with respect to the <j> variable but will be confined to an arbitrarily small window around 6>o 
in the 6 variable. Subsequent applications of G will then give an equal probability of increasing or decreasing the polar angle, 
corresponding to a random walk in 8. 
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FIG. 2: Results of numerical simulations of the discrete time implementation of Grovers algorithm with a noisy oracle. The oracle 
noise is determined by the probability distribution p(e) fSection lll Gil . (P) is the average success rate of 100 trials and N is the 
library size. The curves correspond to calculations with error mag nitudes of e rma = 10~' 5 , 10~' 75 , HT 1 , 10" 1 ' 25 , 10" 1 ' 5 , 10" 1 ' 75 , 
where the labeling goes from left to right. 
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